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Abstract 

The product-wavefunction renormalization group method, which is a 
novel numerical renormalization group scheme proposed recently, is applied 
to one-dimensional quantum spin chains in a magnetic field. We draw the 
zero-temperature magnetization curve of the spin chains, which excellently 
agrees with the exact solution in the whole range of the field. 
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Many interesting problems which have been left unsolved mostly con- 
cerned with strongly correlated systems where many-body effects play es- 
sential roles. For these problems conventional many-body-theoretic tech- 
niques, e.g., the mean-field approximation and the perturbation expansion, 
lose their power; these methods often fail in explaining even qualitative fea- 
tures of the system. As promising substitutes, we have direct numerical 
approaches which are rapidly developing in accordance with the remarkable 
advancement of the computer technology. Examples are the Monte-Carlo 
simulation, the molecular dynamics and the exact diagonalization of Hamil- 
tonian matrices or transfer matrices, which have been successfully applied 
to various problems. 

In these direct numerical methods we face another problem: Since the 
system size N (number of lattice points, for example) which we can deal 
with is always finite, we need extrapolation with respect to A'" to extract 
the behavior in the thermodynamic limit {N — s- oo). As a standard pre- 
scription we have the finite-size scaling method Q, which has been one of 
most successful extrapolation schemes for studying the critical phenomena. 
In applying the finite-size scaling, we must compute the A^-dependence of 
physical quantities for large enough A^, which often requires long compu- 
tation time. For the exact diagonalization method, the large-A^ problem 
is rather serious due to the exponential growth of the Hamiltonian matrix 
dimension, which imposes severe limitation on feasible system size A^. 

The density-matrix renormalization group (DMRG) method invented by 
S. R. White is a very important step to overcome the large- A^ difficulty. The 
method was introduced as an approximate diagonalization scheme and has 
been applied to one-dimensional quantum systems, with results demonstrat- 
ing surprising efficiency of the method 0. The method is also applicable 
to two-dimensional lattice statistical models for diagonalization the transfer 
matrices instead of the Hamiltonian matrices 

Quite recently, two novel algorithms are devised by T. Nishino and K. 
Okunishi, which are closely connected to the DMRG. One is the corner- 
transfer-matrix renormalization group (CTMRG) method |Q, ^ and the 
other is the product- wave-function renormalization group (PWFRG) method |^. 
In the CTMRG, a systematic renormalization scheme for large-lattice-size 
two-dimensional classical spin systems is given in terms of Baxter's corner 
transfer matrix 0]. In the PWFRG, one diagonalizes the "wavefunction 
matrix" instead of the density matrix to obtain the "projection matrix" for 
retained block-state bases and sets up a recursion relation for the "projection 
matrix". As has been shown in refs ^, both methods are highly effi- 
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cient for 2D classical systems. Hence, if we can extend them to ID quantum 
systems, similar efficiency may be expected. For ID quantum systems, the 
PWFRG seems to be more easily implemented than the CTMRG, because 
the latter takes full advantage of the two-dimensionality of the system. The 
actual implementation, however, is non-trivial because the transfer-matrix 
multiplication employed in the original PWFRG lose meaning for quantum 
systems. 

The aim of this article is to show that the PWFRG can actually be ex- 
tended to quantum systems and to demonstrate the efficiency of the "quan- 
tum PWFRG" by applying it to ID quantum antiferromagnets in the uni- 
form magnetic field to obtain the magnetization curve at zero temperature. 
The importance of the study of the magnetization curve itself is clear. First, 
the magnetization process or M(magnetization)-/f (magnetic field) curve is 
a directly and accurately observable quantity by experiments ^ ^. We 
need reliable theoretical calculation on various models for comparison with 
experiments. Second, a magnetized state which is the lowest energy state 
for H > He {He- lower critical field ), is an excited state at H = 0. This 
means that the magnetization curve contains a considerable amount of in- 
formation about the whole energy-level structure of the system. Third, for 
gapless systems, the efficiency (or inefficiency) of the numerical renormal- 
ization groups like the DMRG and the PWFRG has not been fully tested 
yet. Calculations of the magnetization curve would also serve as this test, 
because the finitely-magnetized ground states of Heisenb erg-like quantum 
antiferromagnets are gapless states. 

Consider a spin-5 antiferromagnetic spin chain in a field whose Hamil- 
tonian H is expressed as 

n = J2his„Si+i)-HY,s', (1) 

i i 

where Si is the spin operator at the site i and h is the nearest-neighbor 
coupling function (or local Hamiltonian). We take a unit where g^B = 1 id- 
g-factor, fiB- Bohr magneton) or these factors are absorbed into the field H. 
The "pure" Heisenberg model corresponds to h = h{si,Si+i) = —Jsi ■ Sj+i 
(J: an exchange coupling constant), but we consider more general form of 
h conserving the total S^. Our problem is to find the lowest energy state of 

n. 

Let us briefly explain the PWFRG which is a close relative of the DMRG. 
Recall the infinite-system algorithm of the DMRG ||2|. Write the ground- 
state wave function for 2A^-site system under the free boundary condition, 
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obtained by diagonalizing the Hamiltonian matrix, as 

^W{a\P) (2) 

where a (resp. (3) is the block-state index for the left-half (resp. right-half) 
N sites. For simplicity, we restrict ourselves to the cases where 7i is mirror 
symmetric and the ground-state wavefunction is of even parity with respect 
to the space reflection. If we extend the system by adding two sites at the 
center of the system, the ground-state wavefunction of the two-site extended 
system (with the extended Hamiltonian) can be written as 

^l''+\a,i\j.fi) (3) 

where i and j denote spin states (e.g., i = —S, —S -|- 1, . . . , 5 in the S"^- 
diagonal representation) for added sites. In the DMRG, we form the density 
matrix p as 

p{i', a'\i, a) = ^ ^l''^\a', /3)^'^^+2 .|^.^ 

and diagonalize it to choose m eigenstates with eigenvalues {/j, = 1,2,... , m, 
in descending order in magnitude). We retain the m eigenstates as new block 
states (bases) forming a truncated basis set for (A^ -|- l)-site system. These 
two, (PI) the diagonalization of Tl for the extended system, and (P2) the 
choice of retained bases through the density-matrix diagonalization, are the 
key processes in the DMRG. 

In the PWFRG, these key processes are changed as follows. The process 
(PI) is replaced by, (PI') the improvement of the "input" wavefunction. In 
2D classical systems for which the original PWFRG is developed, the im- 
provement is made by the multiplication of the transfer matrix Q. For ID 
quantum systems, the transfer-matrix multiplication is successfully substi- 



tuted by the modified Lanczos operation ||T^, |ll|] on the wavefunction (see 
eq.(^ below). As for the process (P2), the diagonalization of the density 
matrix is replaced by, (P2') diagonalization of the "wavefunction matrix" 

i>{a,i\p,j) =^Gia,i\j,^), (5) 

where we have suppressed the superscript denoting the number of lattice 
sites. In terms of the "projection matrix" R{a,i\p) which projects the state 
(q, i) into the new block state p, we have 

i>{a,i\/3,j) =Y.R{a,i\p)io^R{(3J\p). (6) 
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Since p = the density- matrix eigenvalues {A^} relate to the wavefunction- 
matrix eigenvalues {w^} as = ^fi @]- Hence the choice of retained bases 
according to the magnitude of is equivalent to that made in the DMRG. 
As in the DMRG, the projection matrix R is also used in the construction 
of the new block Hamiltonian for the size-extended system. 

There is one important process added in the PWFRG: (P3') We set up 
a recursion relation for R as follows. By A{a,i\fi) we denote the "refined" 
projection matrix associated with the improved wavefunction obtained by 
the process (PI')- The recursion relation reads Q 

Rncw{a,i\l^) = Mtj\<^)Ro\d{C,j\r])Mv,i\l^)- (7) 

Updated wavefunction matrix is then given by 

^new(a,i|/3,j) =YRnew{a,i\n)uJf,Rnew{(3,j\fJ-), (8) 

which we put again into the process (PI') so that we can complete one 
iteration sequence. 

We summarize the algorithm of the PWFRG for ID quantum system as 
follows: 

Step 0. Prepare an input wavefunction ^oid and an associated initial projec- 
tion matrix -Roid- Choice of the initial ^'om is somewhat arbitrary; 
the exact ground-state wavefunction for a small-size system is con- 
veniently chosen. Diagonalize the wavefunction matrix to obtain the 
initial projection matrix {Ro\ii{fi\i,a)}. 
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Step 1. Apply the modified Lanczos operation L to ^'oid to obtain the improved 
wave function \I'imp- The operation L = L{a) on a state vector j^* > 
is defined as 

L(a)|^' > = ' ' 



^l + a'^ < x\x > 
= l^'(a) > 
\X> = AH\^ > 

AH = H- < ^[Fl^ > . (9) 
The optimum choice a = a* minimizing < ^(a)\H\^{a) > is given 

by 



«* = (^x3-Vxi + 4xij/(2xi) 

Xn = <^\{AHr\^>. (10) 

Write L{a*) as L*. In most cases we adopt A;- fold operation with 
k>2: 



ll-imp > 



L* 



l^old > 



L*iL*i---iL*\^ old >))■■■) (A: times). (11) 



Step 2. Perform an eigenvalue decomposition of the wavefunction matrix {^'imp(o, i\f3,j)} 
associated with the improved state vector l^'imp >: 

4'iinp(a,i|/?,j) = ^A(a,i|^)w^A(/3,j|/u), (12) 

where the new block index fj, (= 1,2,.. .) is introduced in descending 
order of \u)^\ {{loiI > |c<;2| > •••)• We then truncate the number of 
retained bases to m. 

Step 3. Extend the system size by adding two sites at the center of the chain. 

Then using the "refined" projection matrix {A{a,i\fi)} (/U = 1, 2, . . . , m), 
update the left and right block Hamiltonian matrix just as in the 
DMRG. 

Step 4. Update the projection matrix following eq.(0) and form the new wave- 
function ^'new according to eq.(^). Rename ^new as ^oW, then go to 
Step 1. 
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After repeating the above steps sufficiently many times , we obtain the "fixed 
point" wavefunction. Physical quantities are calculated as the expectation 
value with respect to the fixed-point wavefunction. In particular, the per- 
site magnetization M is calculated as the expectation value of the center 
spins. 

Let us give a comment on the relation between the PWFRG and the 
DMRG. As regards the diagonalization of the wavefunction matrix, a sim- 
ilar process (singular-value decomposition, SVD for short) appears in the 
original DMRG The associated projection matrix R is an important 
object in both methods. Main difference in the two methods lies in the 
updation process of R. In the DMRG, the updation i?oid — > Rnew is "di- 
rectly" given by the SVD of the updated ground-state wavefunction of the 
updated Hamiltonian matrix. Whereas, in the PWFRG, the updation is 
given somewhat "indirectly" through the recursion relation eq.(|^. We do 
not determine -Rncw from the updated wavefunction, but we define Rnew by 
eq.(^ and use it to form the updated wavefunction. The updated wavefunc- 
tion is, then, a good guess of the ground-state wavefunction of the system 
with the extended size (by two sites). We then refine the wavefunction 
by the modified Lanczos operation eq.(^) to obtain the "refined projection 
matrix" A. As the iteration proceeds, the system size becomes larger and 
larger, and at the same time, the wavefunction becomes closer and closer to 
the ground-state wavefunction. The recursion relation eq.(0) implies that, 
at the fixed point (=infinite iterations), R = A should be satisfied, and this 
is precisely the condition for the ground-state wavefunction which must be 
invariant under the modified Lanczos operation. 

We apply the PWFRG to several models of quantum spin chains to draw 
the zero-temperature M—H curves. We present the results in the followings. 



(a) S = 1/2 XY antiferromagnetic chain ||T^, 13| 

The local Hamiltonian h (see eq.(|l])) for this model is given by 

h{s,,sj) = \JMs] + slsJ), (13) 

with Si being |x(Pauli spin matrices). This model is exactly soluble via the 
Jordan- Wigner transformation; the magnetization curve is known to be 

M =-sm-\^). (14) 

IT J 
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This model is a gap less system even aX H = 0, which is reflected in the 
behavior of the M—H curve : 

limM^o+ H{M) = 0. In Fig. | we compare the 
PWFRG calculation with the exact solution, where an excellent agreement 
is seen. It should be noted that even a small number {m = 16) of retained 
bases gives the quite accurate M — H curve. 

(b) Takhtajan-Babujian model 0, ||] {S = 1) 

This model belongs to a family of models called bilinear-biquadratic 
chains whose local Hamiltonians have the form 

h{si,Sj) = \ J\[si ■ Sj + (3{si ■ Sj)^], (15) 

where Si and Sj are S = 1 spin operators. The Takhtajan-Babujian model 
corresponds to /? = —1, which is exactly soluble and is known to be gap- 



less 1 14, 15 1 . The M — H curve can be drawn by solving the Bethe ansatz 
integral equation for the "two-string" (two-down-spin bound state) root den- 
sity to obtain the ground-state energy density with fixed per-site magneti- 
zation M. For general M the solution can only be obtained numerically, 
by converting the integral equation into a matrix equation. As can be seen 
in Fig. H, our PWFRG calculation reproduces the exact curve within a sat- 
isfactory precision even with the small number m (= 20) of the retained 
bases. 

(c) Affleck-Kennedy-Lieb-Tasaki (AKLT) chain |1 (S = 1) 



This model also belongs to the family of the S = 1 bilinear-biquadratic 
chains, whose local Hamiltonian being given by eq.(|T5|) with /3 = 1/3. The 
ground state at = is exactly known to be singlet and gapped (the 
excitation gap A > 0), and is in the "valence bond solid" (VBS) state 
whose wavefunction is the product of matrices with finite dimensions |T7| , p^ . 
The ground-state correlation length and the string order-parameter are also 
exactly known, from which the state is identified to be in the "Haldane 
phase" [|l9|. Increasing the magnetic field H, this VBS state remains to be 
the ground state up to the "lower critical field" Hd above which the ground 
state is magnetized. The critical field relates to the gap A as 

Hci = A, (16) 

which holds for a general class of antiferromagnets so long as A is the singlet- 
to-triplet energy gap. There is another critical field Hc2 (upper critical 
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field) above which the system's per-site magnetization saturates to unity 
(complete ferromagnetic state). For the intermediate region of the field, 
Hci < H < Hc2, the exact ground-state wavefunction has not been known 
even for the AKLT chain. It should be noted that the VBS-state wave- 
function |l^, 18] with finite matrix size m { m = 2 for the AKLT chain 



) is exactly the fixed point of the iterations in the DMRG/PWFRG with 
the number of retained bases just being m [^0|. Hence, by continuity, we 
can expect that the fixed-point wavefunction of the DMRG/PWFRG with 
relatively small m should remains to be a good approximation to the cor- 
rect ground-state wavefunction of the AKLT chain even in the intermediate 
range of the field. 

In Fig. ^ we draw the M — H curve obtained by the PWFRG with a small 
number m (= 20) of the retained bases. The calculation reproduces the exact 
value of Hc2 = 4| J| which can be derived from the stability consideration 
of the saturated state against the one down-spin formation pT| , |2^ , 
Further, the predicted square-root behavior p^, p^, p^] 



M{H) ~ 5 - A2y/Hc2-H ((0 <)A2 : constant) (17) 

is reproduced (Fig. ^). 

As for the behavior near the lower critical field Hd, another square-root 



behavior has been expected p4, M, 26, 



M{H) ~ Aiy^H -Hci ((0 <)Ai : constant), 

Hci = A. (18) 

Although actual studies concerning this behavior have been almost limited 
to the case of the "pure" S = 1 Heisenberg chain without the biquadratic 
term (/3 = in eq.(p!5|)), the theoretical arguments for this behavior are so 
general that we can naturally expect the square-root behavior eq.(^) to 
hold for the AKLT chain either. As can be seen from Fig. |5|, our PWFRG 
calculation is consistent with eq.(|l^); the calculated value of Hd — 0.7| J| is 



in agreement with the previous numerical estimation of A pS, 28]. 



To summarize, in this article we have presented the application of the 
product- wavefunction renormalization-group (PWFRG) method Q to quan- 
tum spin chains to draw the zero-temperature magnetization curve. The 
results show that the method works efficiently even for systems with finitely- 
magnetized ground states which are gapless. 
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We should give a comment on application of the DMRG to the mag- 
netization process of quantum spin chains. In a preliminary stage of the 
present study, we adopted the original infinite-algorithm of the DMRG. We 
then found that the DMRG, applied to quantum antiferromagnets in the 
magnetized region, is very likely to be "trapped" by "metastable states", 
leading typically to oscillatory behavior in the course of the iteration. As 
a result, in most cases, we failed in obtaining the magnetized ground-state 
wavefunction corresponding to the fixed point. In contrast, the PWFRG 
adopted in the present study is fairly stable, always successful in reaching 
the fixed point even in the magnetized region. The above "instability" of 
the DMRG is certainly related to the nature of the present problem. The 
magnetization curve reflects a continuous sequence of level-crossing transi- 
tions induced by varying magnetic field, between the states with different 
values of total . The exact diagonalization (instead of the modified Lanc- 
zos operation used in the present study) in the original DMRG extracts 
the exact lowest energy state of the truncated Hamiltonian at each itera- 
tion step. The state is, however, the lowest energy state of the system at 
that size (~ 2 x iteration steps); owing to the highly degenerate energy level 
structure, there may well occur "level crossings" on the increasing system 
size, leading to the oscillatory behavior in the iteration. 

Of course, we can remove this instability by decomposing the state space 
into subspaces according to quantum numbers to reduce the degeneracy. For 
the problem of the magnetization process the quantum number is the to- 
tal S^{= S^) and we can surely perform the DMRG with fixed S^. This 
approach, however, is rather laborious because what we need is the ground- 
state energy density with fixed magnetization density Sj^ /N (N: lattice size), 
which requires a huge number of fixed-S^ calculations {S^ = 0, 1, 2, ... , NS, 
for the spin S case) for large N, to find the optimum value of Sj^ giving the 
true ground state under the field. On the other hand, with the PWFRG 
method, the system gradually and automatically converges into the "true" 
ground state under the given magnetic field. Due to this "stability", the 
PWFRG will be particularly useful for "critical" systems whose energy spec- 
trum (or transfer-matrix eigenvalue spectrum, for 2D classical systems) is 
highly degenerate. 

The modified Lanczos method used in this article is a two-step restriction 
of the ordinary Lanczos algorithm. Similar restricted-step Lanczos has been 
also employed in a new version of the DMRG (finite-system algorithm) [^] 
which shows very high efficiency. Although precise relation between the 
PWFRG and the new DMRG is not obvious, we can say that the former 
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may, in a sense, be the infinite-system-algorithm version of the latter. 

As for the magnetization process, the "middle- field" phase transitions to 



occur between Hd and Hc2 have been known for some systems pG, 31 



Study of such field-induced ground-state phase transitions is an interesting 
future problem. Finite-temperature behavior |33, 34] of the magnetization 
process is also an interesting problem. 
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Figure 1: M — H curve of the S = 1/2 XY antiferromagnetic chain. 



Figure 2: M — H curve of the Takhtajan-Babujian model. 



Figure 3: M — H curve of the AfHeck-Kennedy-Lieb-Tasaki chain. 



Figure 4: The magnetization curve of the Afileck-Kennedy-Lieb-Tasaki 
chain near the saturation field Hc2- 



Figure 5: The magnetization curve of the Affleck-Kennedy-Lieb-Tasaki 
chain near the lower critical field Hd- 
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